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O ■ 1. INTRODUCTION 

Perturbative Quantum Electrodynamics has been successfully tested to unparalleled 
precision for around half a century. In addition, and rather more recently, the theory's 
behaviour in its strong coupling limit has attracted a growing amount of attention. The 
reason for this is not only that it serves as a simple prototype gauge theory for testing 
nonperturbative calculational techniques which one then goes on to apply to QCD, but 
also because it exhibits interesting behaviour in its own right. In particular, it now appears 
to be quite well established that QED breaks chiral symmetry dynamically as long as the 
coupling is sufficiently large. 

Traditionally, strongly coupled QED has been studied using Schwinger-Dyson, Lattice 
and renormalization group approaches. In this contribution we describe a rather different 
method, namely a variational approach using the worldline formulation of field theory. 
We have developed this technique in the context of a scalar field theory |J and we applied 
it here, for the first time, to a gauge theory with spin-| particles. Although it is too early 
for quantitative results, qualitatively the approach seems promising: The method is for- 
mulated covariantly in Minkowski space-time, the Ward-Takahashi identities are satisfied, 
the connection between the bare- and pole- masses of the electron is gauge independent, 
the known exact result for the infrared dependence of the on-shell renormalized propa- 
gator is reproduced and corrections to the variational approximation can be calculated 
systematically. 

2. THE WORLDLINE TECHNIQUE 

In recent years a rather large literature on the worldline technique has developed |2| , so 
it suffices here to just provide a general outline. We shall concentrate on Green functions 
within quenched QED involving only one fermion line - i.e. the propagator, vertex func- 
tion, Compton scattering etc. The theory is quenched by integrating out the fermions 
and setting the resulting determinant equal to unity. After differentiating the appropri- 
ate generating functional Z with respect to the external fermion's source one obtains the 
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generating functional for the general (2 + n)-point function 
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Here j is the source for the the photon field A, M is the bare mass of the electron, g 
is the coupling constant and Sq[A] is the gauge part of the standard QED action with 
Lagrangian 
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In order to regulate infrared divergences a photon mass m 7 has been introduced and the 
gauge has been fixed covariantly by employing the gauge fixing parameter £. 

We would like to integrate out the photon field. It is not possible to do this at this 
stage because it appears in the electron propagator in Eq. ([]]) in a non-gaussian way. In 
order to proceed, the fermion is re-introduced in terms of a path integral over its worldline 
x(t) as well as a corresponding path integral over a Grassmann valued function ((t) which 
contains the information about the electron's spin degree of freedom || . One obtains 

< V I 7a — I > = dT N ( T ) J d xJ VxV(e iS[x ' (] , (3) 



■up - g4~ M 

where N(T) provides the proper normalization, the integration parameter x is a Grass- 
mannian partner to the proper time T, the four- dimensional bosonic path integral has 
the boundary conditions x(0) = and x(T) = y while the five dimensional fermionic path 
integral has the single boundary condition C M (0) + (^{T) = 27 5 7 M , C 5 (0) + C 5 (^) = ^7 5 - 
The worldline action S[x, (} is given by 

S[x,(] = S [x,(] + 9 I* dt{-x,(t)A»[x(t)} + iF, v [x(t)]C(t)C(t)} ■ (4) 

J 

As expected for a relativistic particle the free part of the action Sbt^C] includes, apart 
from a kinetic term, further contributions which provide coupling between orbital and 
spin motion as well clS db kinetic term for the spin 
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Here m — 0, 1, 2, 3, 5 and the metric g mn = diag(l, —1, —1, —1, —1) is used for short-hand 
notation. Note that the gauge field A appears linearly in Eq. (|j) so that we may now 
integrate out the photons in the generating functional ([!]). Differentiating n times with 
respect to the photon source j (and then setting j = 0) as well as taking the Fourier trans- 
form yields the general Green function of an electron interacting with n external photons. 
With p defined as the ingoing electron momentum and p' ({ft}) the outgoing electron 
(photon) momenta, and the external photon legs as well as the momentum conserving 
delta function removed, one obtains 
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where the T>x implies an integration over the endpoint y in addition to the path integral 
Vx (similarly, we will henceforth define S e fi to be S e s+p'-y). Because the 5 th component of 
the Grassmann function ( is non-interacting, this integration has been carried out in this 
equation. The effective action contains the free term in Eq. (|^) as well as an interaction 
dependent piece 

Si[x,C] = ~f o dhdhj^G^^J^xihlCMJA-Q^ih),^)) 

■ exp [ -iq ■ (x(ti) - x(t 2 )) ] (7) 
Here G^ u is the photon propagator and the fermion current is 

J^[q,x{t),C{t)\ = x„[t) + 2C M (*)g-C(*) , (8) 

which corresponds to the usual Gordon decomposition for the Dirac current of a particle 
of unit mass into a convection and a spin current. 

It is a nice feature of the worldline representation that it is trivial to verify the Ward- 
Takahashi identities for an arbitrary Green function. Because of the Grassmannian nature 
of £(£) a contraction of an external momentum g M into the appropriate current J M just 
yields a total derivative, i.e. % • J, = • x(Tj), which allows one to do the integration 
associated with q i7 i.e. 

[ T dn e**-^) q t ■ J[-q h x(n), (( Ti )} = - (e™ - l) . (9) 
Jo 'i K ' 

The remaining dependence here may be absorbed in a shifted momentum p' in the 

momentum space action S c ff + p' ■ y and so q H ■ G^n i s related to the difference of two 

(2+(n — l))-point functions: 

-q H -G^(j>, v '-{ qi }) = GgtliCp+ay + aKft}) - Gftlx , (10) 
9 

where the notation should be self explanatory. This is the correct result for the Ward- 
Takahashi identity for the (2 + n)-point function untruncated in the external fermion 
legs. Because of the simplicity of the above argument, based only on the structure of the 
fermionic current, it is quite easy to ensure that the variational approximation doesn't 
violate it. 



3. THE VARIATIONAL APPROXIMATION 

Eq. (^j) provides an exact representation of all Green functions occurring in quenched 
QED involving only one fermion line. It is not possible to carry out the remaining path 
integrals over x(t) and ((t) because they do not occur quadratically in the action (see 
Eq. (0)). We shall calculate the Green functions of the theory in an approximate manner 
through the use of the Feynman- Jensen variational principle. This technique has a rather 
long history, having its origin in the polaron problem j| and it was applied, more recently, 
to a scalar relativistic field theory by us . In order to illustrate the idea, let us consider 
the electron propagator (i. e. n = in Eq. (^)). We rewrite the path integral by adding 
and subtracting to S c s a trial action S t : 

J VxV(e^ = (/ VxV( c <&) - (/ VxV( e^) (e*^ , (11) 
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i.e. we need to evaluate the average of the exponential of iAS, weighted with the func- 
tion e lSt . The approximation consists of replacing the average of the exponential by the 
exponential of the average, i.e. 

<e**) A « e'^* , (12) 

where ~ indicates equality at the stationary point of the right hand side under unrestricted 
variations of the functional S t . If one restricts St to depend only on (non-local) quadratic 
functions of x(t) and ((t) then the path integrals may be performed analytically. However, 
only the parameters in the trial action are then left for variation and so it is important to 
note that one can systematically calculate the corrections to the variational approximation 
by including higher terms in the cumulant expansion of which Eq. is the first term. 
Of course, in order for the approximation to be a useful one, the essential physics should 
already be contained in the leading term. 

4. THE VARIATIONAL PROPAGATOR 

In order to assess the feasibility of the program outlined above we have calculated the 
electron propagator, using a very simple trial action 

S t = X'p-y + S [x,C] + [ T dt 1 dt 2 f(t 1 -t 2 ) \ x(h) - x(t 2 )] 2 , (13) 

i.e. one which has the same dependence on the fermionic degrees of freedom ((t) as the 
free action. Only the dependence on the bosonic degrees of freedom x(t) has been changed 
from the free action through the use of the variational parameter A' and the variational 
function / which describes the retardation effects on the electron due to multiple photon 
emission and re-absorption. This is certainly not an optimal trial action and it only serves 
to illustrate some of the qualitative features of what can be expected. In particular, the 
trial action ( |T3"D does not respect the supersymmetry between orbital and spin degrees of 
freedom which is present in the exact wordline action || . Work with a more general trial 
action is underway. Also, we shall only calculate the propagator near its mass-shell: This 
has the great advantage that in Fourier space the trial action effectively diagonalizes, i.e. 

/ dhdtzfih-h) [xih) - x(t 2 )} = ]T A kk ,b k -b' k Y, A kkb\ ,(14) 

J ° k,k'=0 fc=0 

where b k are the Fourier components of the path x(t) and A kk become the variational 
parameters taking the place of the function f(ti — t 2 ) (not to be confused with the 
photon field which has been integrated out). 

With this trial action, and employing Pauli-Villars regularization as well as on-shell 
renormalization, the connection between the physical mass M p h ys and the bare mass M 
becomes 

«k {l - (1 - A) 2 - 2 mX \;J [A - X] } - Ml = (15, 

where A = A'/v4 o ( first discussed by Mano || ). The variational parameters should be 
fixed through the variational equation 

* { ^p hys [l " (1 " A) 2 ] - 2 n[A, A] + V[A, A] } = . (16) 
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Here V[A, A] is a functional which essentially arises from the required averaging of Sj (see 
Eq. flTT|)) while Q[A, A] gets contributions from the averaging of So — S t and the quadratic 
fluctuations of the trial action (apart from the Minkowski formulation it is identical with 
the expression given in Ref. |L|). Both V and Q are gauge independent, as they should 
be, and hence the variational parameters are also gauge independent. If, for M = , 
Eqs. ([TJ|) and ( |T6"D have a common solution with M p h ys 7^ then the theory exhibits 
dynamical mass generation. At present it is not known whether the above trial action 
admits a solution such as this. 

The contributions to the residue, on the other hand, may conveniently be broken up 
into a term F* which is proportional to x an d two terms, F q and F* q , which are not. 
(The subscripts g and qq on these functions indicate whether they arise from the g^ v or 
q^q v terms in the photon propagator.) We obtain 
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where N (A) is an overall factor of no interest here. The important thing to note is 
that this expression only contains infrared divergences in the functions occurring in the 
exponential factor, and these may in fact be extracted analytically. Indeed, the only gauge 
dependent function F qq has the closed form 
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resulting (together with the additional infrared divergence contained in Fg) in a residue 
which contains all its infrared and gauge dependent behaviour in the overall prefactor 
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This is the well-known Bloch-Nordsieck result [[T]. 

In conclusion, we have summarized here the preliminary results which have been ob- 
tained in the first application of the polaron variational technique to QED. We believe 
that the qualitative success shown so far warrants a more quantitative study using a more 
realistic trial action. Work on this is underway and will be published elsewhere. 
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